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Рассматривается способ восстановления фазы сигнала, основанный на численном алгоритме, реализующем операцию «сжатие- 
растяжение», используемый при выделении биоиндикационной информации. Приведены результаты численного исследования 
предлагаемого алгоритма. 


Введение 

Биоиндикационное исследование часто основы- 
вается на анализе слоистых или кольцевых структур 
биообъектов [1]. Информация, содержащаяся в кар- 
тине чередующихся полос разной интенсивности, 
может быть использована для восстановления связи 
между биологическим объектом и окружающей сре- 
дой. Восстановление закодированной подобным об- 
разом информации предлагается осуществлять пу- 
тем вычисления фазы сигнала, представляемого в 
виде колебательного процесса [2]. В результате реа- 
лизации данного подхода получают информацию об 
особенностях изменений локальной структуры ин- 
дикатора. Точность восстановления информации за- 
висит от уровня и вида составляющих шума, содер- 
жащихся в смеси с анализируемым сигналом. При- 
меняя процедуры фильтрации, как во временной, 
так и в частотной областях добиваются требуемой 
точности вычисления фазы. В то же время, для сиг- 
налов с широким спектром частот ошибка восстано- 
вления фазы связана с отсутствием аналитичности. 


Введение аналитического сигнала (АС), позволяет 
однозначным образом определить амплитуду, фазу и 
частоту любой вещественной функции времени 
[3-6]. Наиболее важное свойство аналитического 
сигнала - это причинность его спектра. Спектр АС 
располагается по одну сторону от начала координат 
оси частот [7], т.е. АС, соответствующий действи- 
тельной функции, получается путем обнуления од- 
ной половины спектра частот. Эта операция непро- 
тиворечива только в том случае, когда сигнал сфор- 
мирован как двухполосный, как сумма сопряженных 
комплексных функций с непересекающимися спек- 
тральными полосами. Эти комплексные функции 
будут по определению сопряженными аналитиче- 
скими сигналами. Качественный признак наличия 
двухполосности - это большое число мало меняю- 
щихся периодов у колебания, но при наличии шума 
трудно сделать вывод о выполнении данного требо- 
вания. 

Решение данной проблемы возможно с исполь- 
зованием операции «сжатия-растяжения» сигнала во 
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временной области совместно с фильтрацией в ча- 
стотной области. Преобразование «сжатие», обозна- 
ченное символом С, сжимает периоды колебания, 
которые больше некоторого среднего периода, и ра- 
стягивает те из них, которые меньше, а обратное пре- 
образование - «растяжение» О 1 возвращает сигнал в 
исходное состояние. Применение операции «сжа- 
тие-растяжение» позволяет группировать частотный 
спектр сигнала вблизи несущей частоты, рис. 1. 


•р\ 



Рис. 1. Модуль частотного спектра сигнала: а) до, б) после 
осуществления операции «сжатие-растяжение» 

Численная реализация операции «сжатие-ра- 
стяжение» основана на алгоритме сплайн-интер- 
поляции с использованием априорной информа- 
ции о восстанавливаемой фазе. Предположив, что 
функция фазы ср монотонно возрастающая, полу- 
чим, что обратная ей функция <р _1 будет однознач- 
на, и не будет иметь разрывов, что особенно важно 
для ее численной реализации. 

<р(х) = / с %(х), х(і) = (р-\/ с і), 

где /' - частота несущая, х,і - шкала отсчетов пря- 
мая, обратная. 

Введем преобразование «сжатие-растяжение» 
вдоль оси х следующим образом. Пусть Н есть опе- 
ратор последовательного выполнения полосовой 
фильтрации на несущей частоте / с и преобразования 
Гильберта, тогда для сигнала, представленного в ви- 
де колебательного процесса Т/(х)=[1+ц(х)]со8(р(х), 
где 1 +ц(х) - амплитуда сигнала, алгоритм реализа- 
ции операции «сжатие-растяжение» совместно с 
полосовой фильтрацией может быть представлен в 
следующем виде: 

СР + Ц(х)] С08 ф(х) = {1 + /л[(р \/ с 2)] } X 

X С08 Мф 1 (/ 2)] } = С08 /;, 2 + ц( 2) С03 / с 2, ( 1 ) 

Н[С08 / с 2 + Ц(2)С08 />] = 8ІП / с 2, (2) 

С ■’ 8ІП / С 2 = 8ІП <р(х). (3) 


предположить, что величины щ(х), определяемой в 
соответствии с выражением (4), будет достаточно 
для осуществления начального сжатия спектра 


Ф 0 (х) = агсі8 


Ш/(х) | 

Щх) \ 


( 4 ) 


Затем процесс восстановления фазы реализует- 
ся по формуле (5) [2]. 


<Р 0 (х) = агсіё 


н Щср;'(/ С 2)\ 


Ш' (/Л 

7 -$пМ_ 

/с 1 


( 5 ) 


где «=0,1,2,... 

В результате многократного повторения про- 
цесса восстановления находят последовательность 
значений фазы. Эта последовательность будет схо- 
диться к точному решению в том случае, если при 
неограниченном возрастании числа итераций бу- 
дет существовать предел этой последовательности. 
Наличие сходимости определено в ходе отдельного 
численного эксперимента, по результатам которо- 
го при уровне аддитивного шума о= 0,2 для дости- 
жения заданного уровня допустимой погрешности 
(г=0,0001) было произведено от 5 до 9 итераций. 

Для определения эффективности и целесооб- 
разности применения предложенного алгоритма, 
реализующего операцию «сжатие-растяжение», а 
также для определения свойств данной операции 
проведен замкнутый численный эксперимент. В 
рамках эксперимента было осуществлено числен- 
ное исследование качества оценок фазы. Результа- 
ты работы операции «сжатие-растяжение» получе- 
ны при разных уровнях и видах вносимого в исход- 
ный сигнал шума, а также при исходных сигналах с 
разной шириной спектра частот. 


Описание этапов численного эксперимента 


1. Формирование исходного сигнала и сигнала с 
шумом. 

Исходный сигнал задавался в соответствии с 
математической моделью 1/(х)=[1+ц(х)]со8^(х). В 
качестве фазы в модельном сигнале использова- 
лась функция следующего вида: 


<Р„(х) = 


2л:Т(х-Г) 

N 


+Д 81 П 


2лТ(х-\У 

7Ѵ 2 


-Ъ 8ІП 


2я:(х— 1) 
N 


( 6 ) 


Как следует из уравнений (1-3), преобразование 
«сжатие-растяжение» и полосовая фильтрация по- 
зволяют выделить из колебания с монотонной фазой 
гармоническое колебание с частотой / с , для которого 
находится Гильберт-трансформанта, которая затем 
сжимается и растягивается в обратном порядке. 

Возникает вопрос о реализации описанных 
операций. Для того, чтобы определить ср(х), 
необходимо задать некоторую начальную фазу для 
осуществления преобразований (1-3). Можно 


где <р е - начальная фаза, IV- количество отсчётов, 
а, Ъ - постоянные коэффициенты, Т- количество 
периодов колебания, 7=10. 

Шум формировался в частотной области. За 
центр полосы частот шума была принята несущая 
частота исходного сигнала / с . Частотный диапазон 
шума заполнялся псевдослучайными числами, 
имеющими равномерное распределение. Сигнал- 
шум, вводился в исходный сигнал аддитивно или 
мультипликативно во временной области. 
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2. Восстановление фазы из смеси сигнала с шу- 
мом на основе операции «сжатие -растяжение» и 
определение ошибки восстановления фазы. 

Вычисление фазы осуществлялось в соответ- 
ствии с выражением: 

Ѵ(х) 

ф = агсІ§ — — - ± птг, 

Щх) 

где Ѵ(х) - мнимая часть, Ѵ(х)=Щ ІІ(х)\, Щх) - дей- 
ствительная часть сигнала. 

Восстановление функции фазы производилось 
путем последовательного сшивания элементов 
дискретной фазы. Оценка нормированной сред- 
неквадратической ошибки восстановления фазы 
вычислялась по формуле: 

^(А(Рі-А(р) 2 

р - _і=і 

Ф . N _ ’ 

|5>-Ѵ) 2 

где А(р - разность между исходной и восстановленной 
фазами, Аср - среднеарифметическое значение разно- 
сти фаз, ср‘ - разность между исходной фазой и пря- 
мой, проведенной через начальное и конечное значе- 
ние исходной фазы, (р‘ - среднеарифметическое зна- 
чение разности между исходной фазой и прямой. 

Выборочный ансамбль состоял для всех экспе- 
риментов из 100 разных реализаций сигналов для 
каждого уровня шума. Отдельные контрольные 
эксперименты с большим объемом выборки пока- 
зали, что среднее и дисперсия, полученных оценок 
остаются неизменными, что было расценено как 
наличие статистической устойчивости. При дан- 
ном объеме выборки в силу центральной предель- 
ной теоремы среднее арифметическое значение 
ошибки < в ,,> будет распределено по закону, близ- 
кому к нормальному. При проведении статистиче- 
ского эксперимента представляют интерес случаи с 
малым среднеквадратичным отклонением а ер что 
означает устойчивость алгоритма оценивания фазы 
к входным параметрам, вследствие этого среднее 
арифметическое значение как оценка ошибки фа- 
зы е ѵ будет близко к оценке максимального правдо- 
подобия. 

Шум в частотной области занимал полосу от ну- 
ля до удвоенного значения несущей частоты исход- 
ного сигнала / с . В качестве начальной информации 
для осуществления растяжения и сжатия сигнала 
использовалась исходная фаза. Оценка нормиро- 
ванной среднеквадратической ошибки восстано- 
вления фазы рассчитывалась при изменении отно- 
шения шума к сигналу от 0, 1 до 0,7. 

Из результатов эксперимента, представленных на 
рис. 2, видно, что с увеличением уровня шума растет 
среднее значение нормированной среднеквадрати- 
ческой ошибки восстановления фазы и увеличивает- 
ся ее среднеквадратическое отклонение. Сравнение 
оценок нормированной среднеквадратической 
ошибки восстановления фазы при мультипликатив- 


ном (рис. 2, а) и аддитивном (рис. 2, б) шуме позво- 
ляет сделать вывод о том, что наибольшую ошибку в 
восстановление фазы вносит аддитивный шум. По- 
лученные оценки определяют потенциальную точ- 
ность операции «сжатие -растяжение». 
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Рис. 2. Оценки нормированной среднеквадратической 
ошибки восстановления фазы после операции «сжа- 
тие-растяжение» с исходной фазой при шуме: 
а) мультипликативном, б) аддитивном 

При проведении анализа реальных сигналов ис- 
ходная фаза является неизвестной, поэтому для 
осуществления операции «сжатие-растяжение» 
применяют функцию фазы, восстановленную из 
исследуемого сигнала с шумом. 

По результатам эксперимента с использовани- 
ем в качестве исходной информации для осущест- 
вления операции «сжатие-растяжение» восстано- 
вленной фазы (рис. 3, б, для мультипликативного и 
на рис. 3, г, для аддитивного шума) можно сделать 
вывод о том, что исходная информация о фазе ока- 
зывает существенное влияние на ошибку ее восста- 
новления. Ошибка при сравнении с результатами 
первого эксперимента (рис. 2) выше в 3...5 раз. 

Проведено сравнение полученных оценок нор- 
мированной среднеквадратической ошибки вос- 
становления фазы после применения операции 
«сжатие-растяжение», включающей полосовую 
фильтрацию (рис. 3, б, г), и после полосовой 
фильтрации спектра аналогичных сигналов без 
применения операции «сжатие-растяжение» 
(рис. 3, а, в). Данные результаты отражают то, что 
оценка нормированной среднеквадратической 
ошибки восстановления фазы после применения 
операции «сжатие-растяжение» с исходной фазой в 
несколько раз ниже (для мультипликативного шу- 
ма в 3...7 раз, для аддитивного шума в 2...3 раза) по 
сравнению с результатами полосовой фильтрации 
без применения операции «сжатие-растяжение». 

В результате замены исходной фазы на восста- 
новленную в операции «сжатие-растяжение» ее 
эффективность по отношению к полосовой фильт- 
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Рис. 3 . Оценки нормированной среднеквадратической ошибки восстановления фазы при мультипликативном шуме (а, б) и 
при аддитивном шуме (в, г) после применения: а, в) полосовой фильтрации; б, г ) операции «сжатие-растяжение» с 
восстановленной фазой 


рации без применения операции «сжатие-растяже- 
ние» снижена для мультипликативного шума до 
1,5 раз, для аддитивного шума до 1,3 раз. Данное 
исследование позволяет сделать вывод о том, что 
ошибка восстановления фазы с применением опе- 
рации «сжатие-растяжение» будет тем меньше, чем 
точнее априорная информации об исходной фазе. 

Наиболее значимым исследованием является 
определение влияния величины частотного интервала 
шума Р на работоспособность операции «сжатие-ра- 
стяжение». Данное исследование проводилось при по- 
стоянном уровне шума (о=0,4). Частотный интервал 
шума ^варьировался от 0 до 80 при собственной часто- 
те исходного сигнала /=10. Результаты данного чи- 
сленного эксперимента, рис. 4, показывают, что вели- 
чина частотного интервала шума не оказывает влия- 
ния на изменение ошибки при использовании опера- 
ции «сжатие -растяжение» в отличие от полосовой 
фильтрации без применения операции «сжатие -растя- 
жение», при которой оценка нормированной среднек- 
вадратической ошибки восстановления фазы растет с 
увеличением частотного интервала а дд итивного шума. 


Вторым значимым исследованием является при- 
менение операции «сжатие -растяжение» для сигна- 
лов с «широким» частотным спектром, т.к. они наи- 
более соответствуют реальным сигналам, получае- 
мым при анализе биоиндикационных картин [3]. 
Эксперимент проводился для аддитивного шума по- 
стоянного уровня (а=0,4). Варьируемым параме- 
тром в данном эксперименте являлся интервал ча- 
стотного спектра исходного сигнала (рассматрива- 
лись сигналы, для которых Р= 6, 8, 10, 12 при /=10). 

Из результатов проведенного эксперимента 
(рис. 5) следует, что изменение интервала частотно- 
го спектра исходного сигнала не оказывает влияния 
на изменение оценки нормированной среднеква- 
дратической ошибки восстановления фазы при ис- 
пользовании операции «сжатие -растяжение» до тех 
пор, пока выполняется условие Р</„ а при полосо- 
вой фильтрации без применения операции «сжа- 
тие-растяжение» наблюдается постоянное ее увели- 
чение при расширении спектра частот сигнала. 

В результате замкнутого численного экспери- 
мента показано, что при использовании алгоритма, 
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Рис. 4. Оценки нормированной среднеквадратической ошибки восстановления фазы при изменении диапазона аддитивного 
шума после применения: а) полосовой фильтрации; б) операции «сжатие- растяжение» с восстановленной фазой 
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Рис. 5. Оценки нормированной среднеквадратической ошибки восстановления фазы для сигналов с разной шириной спектра 
частот при аддитивном шуме после применения: а) полосовой фильтрации; б) операции «сжатие-растяжение» с вос- 
становленной фазой 

реализующего операцию «сжатие-растяжение» 2. Величина частотного интервала исходного сиг- 
совместно с полосовой фильтрацией: нала не влияет на точность восстановления фа- 

1. Ошибка восстановления фазы сигнала ниже зы при выполнении условия Р</ с . 

при мультипликативном шуме в 1,5 раза, при 3. Точность восстановления фазы не зависит от 

аддитивном шуме в 1,3 раз по отношению к величины частотного интервала шума, 

ошибке восстановления фазы сигнала при ис- Работа выполнена при поддержке Российского фонда фун- 

пользовании полосовой фильтрации без приме- даментальных исследований, проект «Обь» № 05-07-98009. 
нения операции «сжатие-растяжение». 
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